Coherent transients in optical lattices 
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We use a method based on optical coherent transients to 
study the vibrational coherence lifetimes of atoms trapped in 
the potential wells of a near-resonant optical lattice in the os- 
cillating regime. The dependence of the positions and widths 
of the vibrational Raman resonances of the coherent transient 
spectra on the intensity, detuning, and geometry of the lattice 
beams is investigated and the results are compared with the- 
oretical predictions. It is shown that the principal source of 
broadening of the vibrational Raman resonances is the anhar- 
monicity of the light-shift potential wells. We also perform 
time-resolved measurements of the vibrational excitation of 
atoms equilibrating in an optical lattice. 



I. INTRODUCTION 

In optical lattices atoms are trapped in a peri- 
odic array of micron-sized potential wells. Atoms con- 
fined in these potential wells exhibit a rich dynamical 
behaviour which can be probed by a number of tech- 
niques, including probe absorption spectroscopy het- 
erodyne fluorescence spectroscopy || , recoil-induced res- 
onances 0, intensity correlation spectroscopy HH, and 
observation of the redistribution of photons between the 
lattice- inducing beams 0. Using these techniques, the 
quantized nature of the motion of the atoms in the po- 
tential wells has been demonstrated, yielding information 
on the vibrational frequencies of the atoms, the motion 
of wavepackets, and on the decay of coherences between 
the bound oscillator states. Recently, a novel method has 
been suggested and experimentally demonstrated which 
effectively probes the coherences set up between differ- 
ent vibrational states of the trapped atoms ||||. This 
method has a number of advantages over other schemes, 
as we shall demonstrate. Most importantly, the possibil- 
ity of selectively exciting one particular dynamical mode 
of the lattice allows one to fit theoretical predictions to 
experimental spectra easily and reliably. 
In this paper we are predominantly concerned with res- 
onances of the coherent transient spectra corresponding 
to stimulated Raman transitions between adjacent vibra- 
tional levels. The mechanisms responsible for the decay 
of the vibrational coherences and the resulting widths of 
the Raman peaks are still a subject of debate, and a com- 
prehensive study is as yet absent from the literature. In 



what follows, we present the results of our experiments 
conducted on a three-dimensional lattice of cesium atoms 
in the oscillating regime in which the oscillatory motion 
of the trapped atoms is only rarely interrupted by in- 
elastic photon scattering events. The converse case of 
optical lattices operating in the jumping regime has been 
treated by Mennerat-Robilliard and co-workers [nj. In 
our experiments the intensity, detuning, and geometry of 
the lattice were varied independently, and the resulting 
variations of the positions and widths of the vibrational 
Raman peaks were measured and compared with theo- 
retical models. 

In another set of experiments, we loaded an optical lat- 
tice with atoms having an initial kinetic temperature 
higher than the steady-state temperature in the lattice 
and made time-resolved measurements of the vibrational 
frequency as the atoms came to equilibrium in the lat- 
tice. As the distribution of atoms amongst the vibra- 
tional states evolved, the measured vibrational frequency 
was observed to change owing to the perturbations to the 
energy spacings of the bound levels caused by the anhar- 
monicity of the potential. 

This paper is organised as follows. In secti on O the exper- 
imental setup is described, and in section [illj we present 
the results of the experiments. Section |Tv| contains a brief 
discussion of the various mechanisms involved in deter- 
mining the frequency widths of the vibrational Raman 
resonances observed in optical lattices and a comparison 
of the experimental data with theoretical models. Fi- 
nally, section ^ summarizes our findings. 
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FIG. 1. Experimental setup for a coherent transient mea- 
surement. 



II. EXPERIMENTAL SETUP AND PROCEDURE 



light, it no longer interacts with the atoms and serves 
as a local oscillator for the heterodyne detection of the 
light scattered by the atoms in the lattice. The signal 
arising from the beating of the lattice light scattered by 
the atoms and the local oscillator is then recorded with a 
fast photodiode. Assuming that the coherently scattered 
light from the lattice decays exponentially with time at 
a rate T co h, the observed signal (expressed as a voltage 
V(t)) is of the form 

V(t) oc y/ Iprobelcoh exp(-T coh t) cos[(Aj 

ump S)t] 

) , (1) 

where I pro be and I co h are the intensities of the probe and 
the coherently scattered light, respectively, and Aj ump 
is the amount by which the probe frequency is suddenly 
switched at t = 0. A Fast-Fourier- Transform (FFT) is 
then performed on the signal to yield its frequency spec- 
trum. Our spectra are shifted by subtracting the fre- 
quency jump Aj ump , so that a Rayleigh peak (for which 
<5 = initially) appears at zero frequency in our spectra. 
In the case of the ti me-res olved equilibration experiment 
described in section |lll Q the lattice is loaded directly 
from the MOT with a temperature of order 60/iK, and 
the atoms are then allowed to equilibrate in the lattice 
for a variable length of time. After this period the vi- 
brational excitation of the trapped atoms is probed and 
read-out as described above. 



In our experiment, we trap and cool cesium atoms in a 
standard six-beam magneto-optical trap (MOT) operat- 
ing on the 6S 1/2 (F = 4) -> 6P 3/2 (F' = 5) transition 
After a Is trapping period, the magnetic field gradient is 
switched off and the atoms are further cooled for a few 
milliseconds in an optical molasses to a temperature of 
around 10/xK. The molasses beams are then switched 
off and the lattice beams switched on. The geometry of 
the lattice beams is the standard tetrahedron configura- 
tion jl2|,[l3| with splitting half-angles 8 X and 9 y , where in 
our experiment 9 y is fixed and X is variable (see Fig. [l]). 
The lattice operates on the same transition as the MOT, 
detuned to the low frequency side by 3r to 17T, where 
1/r is the excited state lifetime (r = 2vr x 5.22MHz for 
cesium) . The intensity of the lattice beams in our exper- 
iment lies between 1.5mWcm~ 2 and 3.8mWcm~ 2 . In ad- 
dition to the four lattice beams, a weak (< O.lmWcm" 2 ) 
probe beam propagates either along the x- or the z-axis. 
The detuning S of the probe laser from the lattice laser 
can be controlled by an acousto-optic modulator (AOM) 
connected to a tunable radio- frequency source. 
Coherent transient spectra are recorded as follows: After 
the atoms have equilibrated in the lattice for a few mil- 
liseconds, the probe beam is switched on for 20 — 1000/is 
at detuning 5 (see also section IV). After that time, at 
t = the detuning is increased by 0.9MHz within less 
than 1/is. After the increase in the detuning of the probe 
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FIG. 2. Typical experimental spectra of a Rayleigh (a) and 
a Raman (b) resonance obtained by the coherent transient 
method. In each case, a single feature of the complete spec- 
trum is selectively excited. The magnified plot inset in (a) 
shows more clearly the response around the frequencies of 
the vibrational Raman resonances. 



III. EXPERIMENTAL RESULTS 

Some aspects of the dynamical behaviour of atoms in 
optical lattices are better understood than others. The 
vibrational frequencies of the atoms, for instance, can be 
rather accurately predicted and modelled using a simple 
harmonic oscillator picture (for potential wells that are 
deep enough for a tight-binding model to be valid) Jl2| . 
Discrepancies with experiments can be accounted for by 
the anharmonicity of the potential, which can be treated 
perturbatively fl4|| . 

On the other hand, the lifetimes of atoms in these oscilla- 
tor states and of the coherences between different states, 
and hence the widths they give rise to in the coherent 
transient spectra, are not yet well understood. Several 
mechanisms have been proposed that could be responsi- 
ble for the observed widths of these vibrational Raman 
peaks [fD5|-[r7| . In particular, the photon scattering rate of 
atoms in a light field turns out to be much larger than the 
observed widths of the Raman peaks, indicating that the 
decay of coherences between vibrational levels is strongly 
suppressed. The mechanisms responsible for this sup- 



pression will be dealt with in more detail in section IV A . 
In order to carry out a comprehensive study, we inde- 
pendently varied the detuning, intensity, and geometry 
(characterized by the two splitting angles) of the lattice 
and determined both the position and the width of the 
Raman peaks in the coherent transient spectrum by fit- 
ting Lorcntzians to the experimentally obtained spectra. 
Typical experimental results are shown in Fig. 0. The 
signal-to-noise ratio even of a single-shot measurement 
is good, so that little averaging is necessary. The signal 
recorded by the photodiode is, in general, a sum of expo- 
nentially damped beat-notes between the light scattered 
by the trapped atoms and the probe light. When the 
probe is initially tuned to the same frequency as the lat- 
tice beams, one sees a narrow Rayleigh peak and two 
peaks resulting from off-resonantly excited vibrational 
coherences. When detuning the probe from the lattice 
beams by an amount equal to the separation between 
a Raman peak and the Rayleigh peak (i.e. the vibra- 
tional frequency of the atoms), one can selectively excite 
the vibrational Raman coherences alone and hence ob- 
tain a spectrum that only reflects one excitation mode 
of the lattice, in contrast to probe absorption spectra, 
where resonances associated with many modes are su- 
perimposed. It is found that the resulting spectrum can 
be well fitted to a Lorentzian, and in what follows we 
will refer to the width (FWHM) of such a Lorentzian 
fit as the width of the Raman peak. A further advan- 



tage of the method of coherent transients is the fact that 
the preparation and the read-out of the vibrational co- 
herence grating are separated in time. The frequency- 
selectivity in the excitation phase therefore depends on 
the frequency bandwidth of the initial probe pulse, whilst 
the resolution of the spectrum only depends on the read- 
out time. By making the read-out time as short as pos- 
sible whilst maintaining a reasonable resolution, time- 
dependent studies of the evolution of atoms in an optical 
lattice can be performed. This is discussed in more detail 
in section 
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As can be seen in Fig. ||, the Rayleigh peaks are of much 
smaller width than the Raman peaks. These widths yield 
information on the rate of transfer of atoms between 
neighbouring wells and are discussed in detail in [p^|,p|. 
In this paper, however, we shall concentrate our atten- 
tion on the vibrational Raman peaks. 
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FIG. 3. Plot of the the vibrational frequencies (a) and the 
Raman widths (b). The potential well-depth was varied by 
changing either the intensity or the detuning of the lattice 
beams, keeping the other parameter constant. The widths in 
(b) are plotted against the vibrational frequencies. For clarity, 
in (b) only two representative points have error bars. 
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A. Varying intensity and detuning 

For a tetrahedral lattice (see Fig. [I]) in which the po- 
larizations of the four lattice beams are perpendicular to 
the planes spanned by each pair of nearly co-propagating 
beams, one can derive simple expressions for the vibra- 
tional frequencies in a harmonic approximation of the 
potential close to a minimum |13,f4 1 : 



and 



Q, z = uv(cos ' 



cost 



(2) 



(3) 



where E r = huj r is the recoil energy with tu r — 2?t x 
2.07kHz for an optical lattice employing cesium atoms 
and light tuned close to the D2 resonance line and Uq is 
the potential well-depth given by 



Un = HA 



4J// S 



4(A/r) 2 + 1 



(4) 



(I is the intensity of a single beam, I sat = l.lmWcm -2 , 
and A = uji aser — ojq is the detuning of the lattice beams). 
These theoretical predictions have been verified to yield 
the correct dependence on A, J, and 8, but the exper- 
imental results usually fall short of the theoretical val- 
ues by 15-25 percent. In our experiment, we varied the 
potential well-depth Uq (ranging from Uq ~ 25QE r to 
Uq ~ 2100£V) both by varying the intensity (between 
/ = 1.5mWcm -2 and / = 3.6mWcm -2 ) and the detun- 
ing (between A = — 3r and A = — 17r) and measured the 
positions of the Raman peaks of the corresponding coher- 
ent transient spectra along the x- and the z- directions. 
The results of these measurements are summarized in 
Fig. H (a) and confirm the dependence of the vibrational 
frequencies on Uq as given by equations || and [| The 
splitting half-angles in this experiment were 9 X = 41 ± 1 
degrees and 6L = 38 ± 1 degrees, and the ratio of the 
slopes of Fig. I) (a) is 1.83 ± 0.05, which is close to the 
theoretically predicted value of 1.67 ± 0.06. 
The results of measurements of the Raman widths for the 
cases of the probe propagating along both the x- and the 
z-directions are shown in Fig. Ill A| (b) where, for clar- 
ity, most of the error bars have been omitted [pj| . For 
the most part, the width of the Raman peaks and the 
corresponding vibrational frequencies are proportional to 
each other, with Af2 « 0.45f2. Only for frequencies above 
« 70kHz do the Raman widths deviate appreciably from 
this simple relationship. 
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FIG. 4. Vibrational frequencies (a) and Raman widths (b) 
for different lattice geometries. The detuning and intensity of 
the lattice beams are constant. In Fig.(b), the linear depen- 
dence of the Raman widths on the vibrational frequencies is 
clearly visible. 



B. Varying the geometry 

In order to investigate the dependence of the position 
and width of the Raman peaks on the geometry of the 
lattice, the intensity and detuning of the lattice beams 
(and hence Uo) were kept fixed fixed at / = 3.6mWcm -2 
and A = — 7r, respectively, and the splitting-angle 9 X 
was varied. This resulted in a change in the lattice con- 
stants in the x- and z-directions. In our setup, 9 X could 
be varied between 15±1 and 40 ±1 degrees, which accord- 
ing to equations [2] and || should result in a large variation 
in £l x ,y, but only a small variation in f2 z . Figure |J (a) 
shows the results of the measurements of Vl x as a function 
of sin 0^. As expected, there is a linear relationship, and 
again the experimentally measured frequencies lie about 
20 percent below the theoretical values. 
The measurements of the Raman widths as a function 



4 



of geometry are summarized in Fig. [| (b), where AQ, X 
is plotted against fl x . As before, there is a linear rela- 
tionship between the two, with Afl x rj 0.53512;. Within 
the experimental errors, this agrees with the empirical 



constant of proportionality found in section III A 
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FIG. 5. Equilibration of atoms in an optical lattice loaded 
directly from a MOT with T ~ 60/iK, and the change with 
time in the temperature (a) and observed vibrational fre- 
quency (b) of the atoms is recorded. 



was recorded in about 1ms. Figure || shows the results of 
the corresponding temperature and vibrational frequency 
measurements. The temperature of the atoms reaches its 
steady-state value governed by the lattice parameters in a 
few milliseconds p0| . At the same time, the measured vi- 
brational frequency of the atoms increases as the atoms 
are cooled and then approaches as steady-state value. 
These results can alternatively be visualized by plotting 
the measured vibrational frequency against temperature 
(Fig. ||). Evidently, the colder the atoms are, the higher 
is the observed vibrational frequency. This can be under- 
stood in terms of the anharmonicity of the lattice poten- 
tial, as will be explained in more detail in section [V B| . 
Briefly, the hotter the atoms arc, the larger the popu- 
lation of higher-lying (and more anharmonically shifted) 
vibrational states in the potential wells will be, which 
results in a lower observed vibrational frequency. 
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FIG. 6. Variation of the observed vibrational frequency of 
the atoms with temperature. The dotted line is a linear fit to 
the data. 



IV. THEORETICAL DISCUSSION 



C. Time-resolved probing of the lattice evolution 



A. Mechanisms contributing to the width of the 
vibrational Raman peaks 



In most of the work previously carried out on the ex- 
citation modes of atoms in optical lattices, the atoms 
were allowed to reach equilibrium in the potential wells 
before measurements were performed. To demonstrate 
that we can use the coherent transient technique for 
time-resolved, high-resolution spectroscopic studies, we 
loaded our lattice directly from the MOT, skipping the 
molasses phase. The temperature of the atoms released 
from the MOT was about 60/iK. The lattice was then 
switched on for a variable amount of time, after which 
either the temperature was measured (using the time- 
of-flight method fll9[]), or a coherent transient spectrum 



As mentioned before, whereas the dependence of the 
vibrational frequencies on the lattice parameters is well 
understood and can be theoretically modelled, the widths 
of the Raman peaks of the spectra have so far resisted a 
comprehensive explanation. In this section, we shall sum- 
marize the models proposed so far and discuss them in 
the light of our findings. Also, we give a brief overview of 
the results obtained by other groups and compare them 
to our own experimental results. 

The first rate to consider when dealing with the life- 
times of vibrational levels in optical lattices and coher- 
ences between those levels is the photon scattering rate 
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r' = TUo/hA. After the first observation of distinct Ra- 
man resonances in probe-absorption spectra it was soon 
realized that the lifetimes could not be determined by T' 
since that would have meant that (for most experimental 
parameters) the resonances would not have been resolv- 
able. The experimentally observed widths, however, were 
more than an order of magnitude smaller than those pre- 
dicted by the photon scattering rate. Courtois et al. |15(| 
proposed the Lamb-Dicke effect |2lJ arising from sub- 
wavelength localization of the atoms as an explanation 
for the unexpected narrowness of the resonances p2[ . 
According to this theory, the rate of inelastic scattering 
events in which the vibrational quantum number changes 
is reduced by the Lamb-Dicke factor (n + l)E r /%Q, in 
favour of clastic scattering events involving no change 
in vibrational level. For typical parameters, the Lamb- 
Dicke factor is of the order of 1/10 or less, so that the 
inelastic scattering rate F" oc T' (n+\)E r /Ml comes close 
to the experimentally observed widths. In experiments 
involving shallow potentials with small Lamb-Dicke pa- 
rameters, however, Kozuma et al. still found lifetimes 
that were much longer than predicted by the Lamb-Dicke 
model. Phillips et al. interpreted this as evidence for the 
existence of the transfer in inelastic scattering processes 
of coherence between different pairs of vibrational lev- 
els. Also, results of other experiments suggested that 
the broadening of the vibrational Raman resonances was 
dominated by the anharmonicity of the lattice potential, 
which was further demonstrated through the observation 
of revivals in wavepacket oscillations. In the following, 
we will show that our experiments support that inter- 
pretation. Before discussing the experimental findings in 
detail, however, we shall give a brief outline of the effects 
of anharmonicity in optical lattices. 



B. Effects of anharmonicity 

When deriving the expressions for the vibrational fre- 
quencies of atoms in optical lattices (equations |^ and |^) , 
the true light-shift potentials are often approximated by 
harmonic potentials characterized by the curvature of the 
lattice potentials close to a local minimum. The effect of 
the anharmonicity of the potential can be treated pertur- 
batively by considering the correction due to the fourth- 
order term in the expansion of the potential in first-order 
perturbation theory. Taking into account the geometrical 
factor sin 9 X characterizing the lattice constant in the re- 
direction, this leads to the vibrational frequency fi n)n +i 
associated with a pair of adjacent vibrational levels char- 
acterized by n and n + 1 being reduced with respect to 
the harmonic frequency f2o so that 



Qn,n+i ~ SIq — u> r (n + 1) sin 2 



(5) 



(as has been shown experimentally the effect of 

the anharmonicity on the experimentally observed Ra- 
man resonances can be estimated. Firstly, the observed 
vibrational frequency can be defined as an average over 
all the anharmonically shifted pairs of vibrational levels, 
weighted by their respective Boltzmann factors, leading 
to the following expression ]T^ |: 



n w o (i - 2C 



(6) 



This indicates that all the vibrational frequencies should 
be reduced by a constant fraction, and that that frac- 
tion will be proportional to the temperature T of the 
atoms in the lattice for a constant well-depth Ua- This 
satisfac torily accounts for our observations presented in 
section 
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Likewise, the broadening of the resonance due to the pres- 
ence of several vibrational frequencies can be modelled by 
assuming that only levels with a Boltzmann factor (rela- 
tive to that for n = 0) larger than some small, arbitrary 
number b contribute appreciably to the observed spec- 
trum. The most shifted vibrational frequency, which is a 
measure of the frequency spread, is then given by 

An^sme x Cln{b- 1 )Lu r \^- = ^-\n{b- 1 )n Q . (7) 
V E r 2 

This is proportional to the vibrational frequency Qq for 
any value of 6, so that we would expect that the true 
linewidth, calculated from a properly weighted sum over 
all vibrational levels, also to be proportional to £Iq. This 
is precisely what we observe in our experiments. Figure |?j 
illustrates the dependence of the vibrational frequency 
and width of the Raman peak on the geometrical factor 
sin 6 X . In a more quantitative treatment, we have per- 
formed numerical simulations of the coherent transient 
spectra taking into account the frequency spread by su- 
perposing the contributions of all possible pairs of vibra- 
tional levels. The results exhibit the same characteristics 
as the qualitative models described above. 
In principle, the anharmonicity of the potential should 
also manifest itself in an asymmetry of the Raman peak. 
One expects a slower drop-off of the Raman peak on the 
low frequency side because the anharmonic shifts are all 
in the same direction (see equation^). This effect, how- 
ever, will only be noticeable if the frequency spread due 
to anharmonicity is appreciable compared to the widths 
of the individual contributions, which are caused by the 
dissipative mechanisms described above. In our experi- 
ments, we have not found any clear evidence of an asym- 
metry of the Raman peaks. In Fig. 7], however, the 
Lorentzian fit to the broader Raman peak (b) has a larger 
% 2 -value then the narrower peak (a), which suggests that 
the former is distorted by anharmonicity. 



Assuming that in an optical lattice the vibrational lev- 
els are populated according to a Boltzmann distribu- 
tion H|, and that the ratio C — ksT '/2Uq is a constant 
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FIG. 7. Example of Raman peaks for two different geome- 
tries of the lattice with identical detuning and intensity of the 
lattice beams. When the splitting-angle is increased, both the 
position and the width of the Raman peak change. The dot- 
ted lines are Lorentzian fits. 



C. Comparison of theory and experiment 

One of the main objectives of our experiments was to 
determine the dominant mechanism involved in deter- 
mining the width of the Raman resonances in the coher- 
ent transient spectra. In particular, we were interested 
in establishing whether the Lamb-Dicke model correctly 
predicted the dependence of the widths on all the experi- 
mental parameters. Since our experiment was performed 
on a 3-D lattice, whereas the Lamb-Dicke treatment of 
Courtois et al. is based on calculations in one dimension, 
we have generalized that derivation to take into account 
the dependence of the lattice constant on the geometry 
of the lattice. This leads to the appearance of an extra 
geometrical factor sin 6 X in the expression for the rate 
r"(n) of inelastic Raman processes leading to a change 
in the external state of the atom along the x-direction 



T"(n) ocT'sm9 x {n + l)J — 

V u 



(8) 



This reproduces the observed linear dependence of Afl o n 
fi when the splitting-angle 6 X is varied (see section III E| ) . 
It also turns out that when the lattice intensity is varied, 
the expression for T" reduces to T" oc f2. When Uq is 
varied by changing the detuning, however, the relation- 
ship between Q, and T" (and hence Af2) is expected to 
be of the form T" oc f2 3 (for large detuning), which is not 
at all what we observe. In fact, upon closer inspection 
the values of A£l for identical f2 but different detuning 
do differ slightly from each other, the value of Ail for 
the larger detuning being smaller, as expected. The dif- 
ference, however, is very small and does not reflect the 
predicted f2 3 -dependence. This is perhaps only to be ex- 
pected because, for the case of the largest detunings used 
in these experiments, the proximity of the neighbouring 
F = 4 — > F' = 4 transition means that our approxi- 
mation of an isolated transition is of dubious validity. 
Nevertheless, we find that the observed dependence of 
Af2 on f2 agrees very well with broadening due to anhar- 
monicity, as derived above. These observations suggest 
that the main broadening mechanism is, indeed, the an- 
harmonicity of the potential, and the dissipative mech- 
anisms based on inelastic photon scattering are strongly 
suppressed. In the harmonic limit, where the homoge- 
neous width is small compared to the anharmonic fre- 
quency difference between successive vibrational levels, 
the suppression of line-broadening because of the trans- 
fer of vibrational coherence ]24j] between pairs of levels 
in spontaneous Raman processes is expected to be ap- 
preciable. The lifetimes of the vibrational coherences are 
then determined by the cooling rate, rather than by the 
inelastic photon scattering rate, which in the oscillating 
regime is given by T coo i — ATE r /TiA |i~7| and, in our ex- 
periments varies between 3kHz for small detuning to less 
than 1kHz for the largest detuning employed. For large 
detuning, the associated width is given by Ail oc f2 2 . 
Conversely, in the anharmonic limit, it is the spatial lo- 
calization that is expected to give rise to the suppres- 
sion of inelastic Raman processes. Clearly, the interplay 
between the transfer of coherence, the Lamb-Dicke nar- 
rowing and the line broadening mechanisms is hard to 
model quantitatively for an intermediate case. Never- 
theless, our results strongly suggest that the broadening 
associated with the anharmonicity generally dominates 
the widths of the Raman peaks, with dissipative mecha- 
nisms contributing at a much reduced level. 



V. CONCLUSION 

We used the method of coherent transients in order to 
measure the vibrational frequencies of atoms in optical 
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lattices and the widths associated with the corresponding 
Raman peaks of the coherent transient spectra. Vary- 
ing the intensity, detuning, and geometry of the lattices 
beam, we confirmed the validity of the theoretical predic- 
tions for the vibrational frequencies based on a harmonic 
oscillator model of the potential wells. Our measure- 
ments of the Raman widths as a function of intensity, 
detuning, and geometry all yielded a linear dependence 
on the vibrational frequencies, which strongly suggests 
that anharmonicity is the main broadening mechanism. 
This implies that the decay of the coherence between ad- 
jacent vibrational levels is suppressed, which we ascribe 
to the combined effects of the transfer of coherence be- 
tween pairs of vibrational levels in inelastic Raman pro- 
cesses and the suppression of the inelastic rate arising 
from the sub-optical wavelength spatial localization. The 
anharmonicity was also found to be responsible for a re- 
duction in the observed vibrational frequencies relative 
to the predictions of a harmonic model, which was tested 
by measuring the vibrational frequencies of ensembles of 
atoms with different temperatures in otherwise identical 
optical lattices. 
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